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ABSTRACT 



We present high angular resolution (down to (X'3 ~ 13 AU in diameter) Submil- 
limeter Array observations of the 880 fim (340 GHz) thermal continuum emission from 
circumstellar dust disks in the nearby HD 98800 and Hen 3-600 multiple star systems. 
In both cases, the dust emission is resolved and localized around one stellar compo- 
nent — the HD 98800 B and Hen 3-600 A spectroscopic binaries — with no evidence for 
circum-system material. Using two-dimensional Monte Carlo radiative transfer calcu- 
lations, we compare the SMA visibilities and broadband spectral energy distributions 
with truncated disk models to empirically locate the inner and outer edges of both disks. 
The HD 98800 B disk appears to be aligned with the spectroscopic binary orbit, is in- 
ternally truncated at a radius of 3.5 AU, and extends to only 10-15 AU from the central 
stars. The Hen 3-600 A disk is slightly larger, with an inner edge at ~1 AU and an outer 
radius of 15-25 AU. These inferred disk structures compare favorably with theoretical 
predictions of their truncation due to tidal interactions with the stellar companions. 



Subject headings: circumstellar matter — planetary systems: protoplanetary disks 
stars: pre-main sequence — binaries: spectroscopic — binaries: visual 



1. Introduction 



A large fracti on of the stars in the solar neighborhood are members of g ravitationally bound 



multiple systems (lAbt h Levy 



1976 



Duquennoy Mayor 



1991 



Lada 



mult i plicity fractions are found in nearby young stellar associations (IGhez et al 



2006). Similar or higher 



1993 



Leinert et al 



1993 



Simon et al.l ll995). confirming that multiplicity is a common outcome in any star formation 



environment. With such a substantial number of stars residing in multiple systems, it is natural 
to also consider the potential impacts of stellar companions on the planet formation process. The 
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structures of the circumstellar disks used to build p lanets can be strongly a f fected by tidal interac- 
tions with the stellar compon ents in these systems (|Lin &: Papaloizoulll993l ; lArtvmowicz Lubow 
1994 ; iPierens fc Nelson! 120081 ), These interactions can open gaps in disks, truncate their edges, and 
potentially expedite their dissipation before planetesimals can be made. But remarkably, despite 
these formidable dynam ical hazards, a relatively large fraction (~25%) of exoplanet host stars 
have stellar companions (jRaghavan et al.ll2006l : iDesidera Barbieril 120071 . and references therein). 
Observations of circumstellar material in young multiple star systems can be used to empirically 
characterize the impact of these tidal interactions on disk structures, and therefore to place con- 
straints on when and where planet formation can potentially be accommodated in the presence of 
a stellar companion. 

The fate of the circumstellar material in a young multiple star system is primarily depen- 
dent on the orbital separatio n (a), eccentricity (e), and mass ratio (yu) of the stellar components 
(jArtymowicz fc Lubowlll994l ) . Systems with large separations (a ~ hundreds of AU) impart little or 
no dynamical effects on their circumstellar material, leaving disks around each stellar component 
that are similar to those around single stars. Conversely, individual disks in a small-separation 
system (a ~ a few AU or less) will likely not survive. Instead, these systems can host a single 
circum-system disk with a dynamically cleared central cavity out to a radius comparable to the 
stellar separation (~2-3a). Qualitatively, such theoretical predictions have been broadly supported 
by observations. The disks around small-separation spectroscopic binaries are like those around 
single stars except for a diminished dus t emission signal shortw ard of ~5^tm, consistent with a 
lack of warm material inside a few AU (| Jensen fc Mathieulll997l ). Likewise, disks in wide binary 



systems are comparable to thei r more isolated counterparts, although the secondary disk signatures 
are o ften weak or absent (e.g., I Jensen fc Akesonl 120031 : 1 Andrews Williams! l2005al : IPatience et al 



20081 ) 



But most multiple systems in b oth the field and you ng clusters have orbits with intermediate 



separations (a ~ tens of AU; e.g., iMathieu et al.l 12000 ) . The disks in these systems suffer the 



most dramatic effects of star-disk interactions, resulting in their external truncation at a fraction 
of the component separation (~0.2-0.5a), or perhaps their complete dispersal. Single-dish radio 
photometry surveys have demonstrated that the ~1 mm luminosities for medium-separation binaries 
are si gnificantly lower than their more closely or widely separated counterparts (as we ll as single 
stars; Jensen et al. 1994 . 1996a : Osterloh <fe Beckwith 19951: Andrews Sz Williams 2005b ). Since the 
dust emission at these wavelengths is optically thin ( Beckwith et al. 199ol ). this luminosity deficit is 
interpreted as indirect evidence for a diminished disk mass due to the tidal truncation of the outer 
disk. However, this total mass deficit does not necessarily have any bearing on the potential for 
making planets around medium-separation multiples. The likelihood for planet formation depends 
more on reaching a critical disk density, and therefore it is more important to identify how that 
total mass is distributed spatially in such systems. Doing so involves locating where the disks in 
multiple star systems are truncated using high angular resolution observations of their millim eter 



continuum emission (e.g.. IJensen et al.lll996bl : lAkeson et al. 



1998 



Guilloteau et al.lll999l . 12008ft . 
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In this article, we present the first set of high angular resolution millimeter observations of the 
circumstellar material in two nearby multiple star systems, HD 98800 and Hen 3-600. Both systems 
are members of the TW Hya associati on and are therefore among the neares t pre-main sequence 
stars, with a distance of roughly 50 pc (jKastner et al.lll997l ; IWebb et al.lll999l ). HD 98800 is a rare 
quadruple star system, consistin g of a pair of spec t roscopic binarie s with ~lyr periods separated 
by / /8 (36 AU) on the sky (e.g.. iTorres et al.lll995l : lTokovininlll999l ). Infrare d measurements have 



shown that only the HD 98800 B system harbors a circumstellar dust disk (jKoerner et al. 



Prato et al. '. 


2001 




(Furlan et al 


. 2007) 



2000 



The coronal X-ray line ratios and very weak Ha emission observed from 



HD 98800 B suggest that the stars are not actively accreting gas from the disk (jKastner et al 



2004 



Dunkin et al.l 119971 ) . Hen 3-600 appears to be a hierarchical triple star system, comprised 



Torres et al 



2003) 
19891 ). 



of a primary spectroscopic binary of uncertain period (jMuzerolle et al.l 12000; 
and a secondary visual companion located l'/5 (66 AU) to the southwest (jde la Reza et al 
The signatures of a circumstellar disk, including spectroscopic and X-ray line ratio evidence for 
accre tion onto the central stars , have been identified around the Hen 3-600 A spe ctroscopic binary 
(e.g., Ijayawardhana et al.lll999l : iMuzerolle et al.l 120001 : iHuenemoerder et al.l 120071 ). 



At their advanced ages of ~10Myr, the nature of the circumstellar material in these low- 
mass pre-main sequence systems is a topic of active debate. While the accretion signatures in the 
Hen 3-600 A disk are supportive of a primordial reservoir, it is unclear whether the HD 98800 B disk 
contains gas or has evolved into a debris-only system (see the detailed discussion on this topic by 
Akeson et al.ll2007l ). Regardless, the proximity of these systems offers unusually detailed views of 
their circumstellar structures. We aim to use the high spatial resolution afforded by that proximity 
to characterize how the disk structures have been shaped by their dynamical interactions with 
the stellar components in each system. We present our observations in §2 and a brief qualitative 
analysis of the data in §3. We discuss the radiative transfer modeling used to determine the disk 
structures in §4, and synthesize the results with theoretical predictions for disk truncation in §5. 



2. Observations and Data Reduction 

The HD 988 00 and He n 3-600 systems were observed with the Submillimeter Array interferom- 



eter (SMA: iHo et al.ll2004l ) at Mauna Kea, Hawaii in the compact (baselines of 16-70 m), extended 
(baselines of 28-226 m) , and very extended (baselines of 68-509 m) configurations on 2008 March 
9, February 21, and April 3, respectively. The SMA double-sideband receivers were tuned with 
a local oscillator (LO) frequency of 340.755 GHz (880 /im). Each sideband provides 24 partially 
overlapping 104 MHz spectral chunks centered ±5 GHz from the LO frequency. The correlator was 
configured to subdivide one chunk in each sideband into 128 spectral channels to sample any emis- 
sion from the CO J=3— 2 transition (345.796 GHz) at a velocity resolution of 0.7 km s . Each of 
the remaining chunks was split into 32 coarser spectral channels to observe continuum emission. 



The observations cycled through the two targets and TW Hya (Andrews et al., in preparation) 
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between integrations on the nearby quasar J1037— 295, with a total cycle time of 18 minutes (5 
minutes per each disk target, 3 minutes on the quasar) in the compact and extended configurations 
and 12 minutes (3 minutes for each disk target and quasar) in the very extended configuration. The 
bright quasar 3C 279 was observed every other cycle to help assess the quality of phase transfer. 
Additional observations of bright quasars (3C 84, 3C 273, 3C 279, and 3C 454.3) and Callisto were 
made for bandpass and absolute flux calibration when the targets were at low elevations (<20°). 
The observing conditions were excellent, with zenith opacities ranging from 0.04-0.08 (corresponding 
to 0.8-1.6 mm of precipitable water vapor) and well-behaved phase variations on timescales longer 
than the calibration cycle. The extended array observations (2008 February 21) were conducted 
during a configuration change, and therefore only utilized 4 of the 8 SMA antennas. 

The data were edited and calibrated with the MIR software package^ The bandpass response 
was calibrated with observations of bright quasars, and broadband continuum channels in each 
sideband were generated by averaging the central 82 MHz in all coarse resolution spectral chunks. 
The visibility amplitude scale was set by bootstrapping quasar flux densities from the observations 
of Callisto. The absolute flux calibration uncertainty from this technique is estimated to be roughly 
10%. The antenna-based complex gain response of the system as a function of time and elevation 
was determined with reference to the observations of J1037— 295 (the nearest sufficiently bright 
calibrator), located midway between and ~10° from each of the science targets. The phase transfer 
between the calibrator and science targets over that angular separation is not ideal, especially given 
the low observing elevations for these southern systems. Atmospheric decoherence on these scales 
will add phase noise that acts to smear out the emission from the targets. Typically, observations 
of a second quasar are utilized to estimate the level of decoherence in terms of the size of the 
"seeing" disk. Unfortunately, the next nearest bright quasar, 3C 279, probes the atmosphere at 
a much larger 30-40° separation from the targets. Nevertheless, the observations of 3C 279 were 
used to estimate the quality of phase transfer in an extreme case and showed a seeing disk size 
of ~0 / .'3. The actual phase decoherence between the targets and J1037— 295 should be much less 
severe, producing significantly smaller (but still uncertain) seeing disks. 

After the gain calibration, the independent sets of visibilities from each configuration and side- 
band were compared to check for consistency on overlapping spatial scales. After finding excellent 
agreement, all of the data for each target were combined. The standard tasks of Fourier inverting 
the visibilities, deconvolution with the CLEAN algorithm, and restoration with a synthesized beam 
were conducted with the MIRIAD package. Maps of the continuum and CO J=3— 2 line emission 
were made by naturally weighting the visibilities. The resulting continuum maps have 1 a rms 
noise levels of 3.5 mJy beam -1 , with synthesized beam dimensions of 0792 x / .'68 at a position 
angle (PA, measured east of north) of 11° for HD 98800 and l'.'ll x 0"74 at PA = 18° for Hen 3-600. 
To maximize the sensitivity to line emission, only data from the compact configuration were used 
to generate CO J=3— 2 spectral images, with synthesized beam dimensions of 3'.'1 x l'.'8 at PA = 



See http : //cf a- www. harvard. edu/$\sim$cqi/mir cook. html 



- 5 - 



164° and 4'.'0 x l'.'8 at PA = 168° for HD 98800 and Hen 3-600, respectively. No line emission was 
detected for either target, with 3cr upper limits of ~0.6 Jy beam -1 in a 0.7 km s -1 channel. 



3. Results 



The synthesized images of the 880 fj,m (340 GHz) dust continuum emission in the HD 98000 
and Hen 3-600 systems are shown in Figure [TJ The HD 98800 map exhibits emission with an 
integrated flux density of 110 mJy and peak of 82 ± 3mJy beam" _1 , centered at a = ll h 22 m 05"24, 
5 = — 24°46 / 39 / . / 10 (J2000). Based on the reference coo rdinates and proper m otion estimates from 



Hipparcos (which resolved the HD 98800 A-B system; iPerrvman et al.lll997l ). the SMA emission 
centroid is coincident with the HD 98800 B spectroscopic binary position to within ~50 mas, well 
below the SMA astrometric uncertainty of ~0'. / l (and roughly on the scale of the Ba-Bb separation). 
The location of the HD 98800 A spectroscopic binary is marked on Figure[T]with a star symbol / .'8 to 
the south of the map center. The Hen 3-600 system shows emission with an integrated flux density of 
75mJy and peak of 57 ± 3mJy beam" 1 , centered at a = ll h 10 " 1 27!84, 6 = -37° 31 / 51?71 (J2000). 
Using the Two Micro n All Sky Su r vey ( 2MASS) astrometry (jCutri et al.l 120031 ) and the proper 
motion determined by I Torres et al.l (120061 ) , we find that the SMA emission centroid is coincident 
with the Hen 3-600 A spectroscopic binary position to within 15 mas in right ascension and 0'. / 15 
in declination (the measured position is offset to the north, away from Hen 3-600 B). Given the 
astrometric errors that may occur for such a low declination target as well as the potential confusion 
in the 2MASS catalog astrometry from the close companion, we suggest that the continuum emission 
measured with the SMA is consistent with being centered on Hen 3-600 A. The location of the Hen 3- 
600 B companion is marked in FigureQ]with a star symbol ~1'.'5 to the southwest of the map center. 

An examination of the SMA visibilities confirms that the emission from each target is partially 
resolved (see §4). An elliptical Gaussian fit to the HD 98800 B visibilities indicates an elongated 
emission morphology, with a high aspect ratio corresponding to an inclination i ~ 55-71° at a major 
axis orientation PA ~ 130-159° east of north. This inferred disk geometry is in excellent agreement 
with the much more precisely determined orbital geometry of the HD 98800 B spectroscopic binary, 
for which iBoden et al.l (|2005l ) have measured i = 67° and f2 = 338° (corresponding to PA = 158°). 
Within the uncertainties of the SMA data, the HD 98800 B disk and spectroscopic binary orbital 
plane are aligned. A similar fit to the Hen 3-600 A visibilities suggests a more circular emission 
morphology, with inclination estimates of i ~ 25-45° and a major axis orientation PA ~ 145-185° 
east of north. Lacking orbital information for either the spectroscopic binary or the wider visual 
pair, the relative orientation of the star and disk axes in this system remains unclear. 
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Disk Models 



The 880 fim. visibilities presented here provide the first spatially resolved information about 
the HD 98800 B and Hen 3-600 A circumstellar disks. We aim to use these data to directly locate 
the outer disk edges, where the circumstellar material could be truncate d by tidal interactions wit h 
their nearby stellar companions. Moreover, following the analysis of I Jensen &: Mathieul (|1997l ). 
we seek to indirectly locate the inner disk edges sculpted by their central sp ectr oscopic bina- 



ries b ased on the high-quality infrared spectra presented by lUchida et al.l (J200J) and iFurlan et al 



(|2007l ). To determine these disk edges, we utilize a parametric tru ncated disk structure model a nd 
a two-dimensional Monte Carlo radiative transfer code (RADMC; see iDullemond &: Dominikll2004l ) to 
compute synthetic 880 /im visibilities and spectral energy distributions (SEDs) that can be com- 
pared with observati ons. A detailed description of the modeling methodology was provided by 
Andrews et al.l (|2009l ). and so here we focus only on the key modifications made to accommodate 
truncated disk structures. 

We parameterize the axisymmetric density structure for the disk models based on two different 
radial zones, including a standard "outer" disk and a central disk "cavity" that contains substan- 
tially less material. The surface densities have a £ oc 1/R behavior for i? C av < R < -Rout in the 
outer zone, and a constant value (depleted by a factor <5s from the value at i? C av) inside the cavity 
{Rin < R < R C av)- The vertical density structure is assumed to be Gaussian and is parameterized 
by a scale-height distribution H oc R l+ ^ that can be artificially "puffed" -up by a factor 5h at 
the cavity edge (R CSC v) to account for the local ly increased energy input from the normal ir radia- 
tion of the optically thick outer disk zone (e.g., IDullemond et al.ll200ll ; iD'Alessio et al.ll2005l ). The 
densities in these structures are normalized with a total dust mass, Md us t. 



We use the same set of dust opacities adopted by I Andrews et al.l (|2009h . computed for a power- 
law size distribution up to a maximum grain radius of 1mm. For a small fraction, / sg , of the dust 
mass inside the disk cavities, the maximum grain size was decreased to 1 fim (hereafter "small" 
grains) to more faithfully reproduce the solid-state dust emission features for each target. In these 
models, we assume the only he ating source was irrad iation from a central stellar photosphere, repre- 
sented by appropriately scaled iLejeune et al.l (Il997l) spectral templates that match the component- 



resolved optical and near infrared photometry (|Gregqrio-Hetem et al 



Webb et al 



1999; 



Low et al 



Geoffray Sz Monin 



2001 



ired p 
1999; 



Weintraub et al 



2000; 



Koerner et al 



1992 : iSoderblom et al 
2OO0l : 



Prato et al 



1998; 
2001; 



. Table Q] lists the stellar parameters adopted in the modeling. The corre- 
sponding photosphere templates are also used to help construct the complete broadband SEDs for 
each disk. Since neither the HD 98800 A or Hen 3-600 B components are found to harbor any disk 
material, we subtra c t their stellar temp lates from the composite Spitzer IRS spectra for each system 



(Uchida et al. 



2004 



Furlan et al.ll2007l ) to derive the star+disk spectra for HD 98800 B and Hen 3- 



600 A alone. Assuming that the emission at longer wavelengths is generated solely by the disks, 
we include the fa r-infrared and millimeter photometry from the l i terature in the S EDs without any 
such corrections (jjensen et al.lll996al : ISylvester et al.lll996l . l200ll : lLow et al.ll2005l ). 
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In this modeling effort, no attempt was made to account for irradiation from individual com- 
ponents of the eccentric spectroscopic bin aries or from any external disk irradiation by the close 
stellar companions (e.g.. lNagel et al.ll2009l ). For convenience, th e updated HD 98800 Hipparcos dis- 
tance of 45 pc was adopted for both systems (jvan Leeuwenll2007l ). Note that this kind of parametric 
modeli ng is still subject to some intrinsic degeneracies, which have been discussed at length else - 



where (jThamm et al.1 1 19941 : IChiang et al. 



2001 



Andrews 8z Williams 



2007a 



Andrews et al J 12003 ). 



Some degeneracies specific to this truncated disk structure model are highlighted below. Because 
the 880 fim emission from the target disks is relatively faint, the exploration of parameter-space 
was limited. Focus was placed on estimating R CSLV and R out for reasonable assumptions about the 
ranges of other input parameters. 

The best estimates for the HD 98800 B and Hen 3-600 A disk structure parameters described 
above are compiled in Table [2J The observed broadband SEDs and Spitzer IRS spectra from the 
literature are shown in Figure with the composite stellar photosphere spectra overlaid in gray. 
The best-fit synthetic SEDs for disk models with several fixed values of the outer radius (i? ou t) 
are plotted for comparison. Neither disk exhibits dust emission in excess o f the photosphere at 



1999: 


Uchida et al. 


2004; 


Fur lan et al. 


2007) 



sion peak near 20 /im were reproduced in the models by varying the disk cavity parameters, {Ri n , 
Raw, <5ei 8h} (and to some extent / sg ). Of these, we perhaps have the best handle on i? cav due 
to the prominent role it has in shaping the rising mid-infrared spectrum and setting the wave- 
length of the secondary emission peak. Our estimate of -R C av = 1 AU for the Hen 3-600 A disk is 



in good agreement with the original measurement from these data by lUchida et al.1 (120041 ) . but the 
Rc»v = 3.5 AU infer red here for the HD 98800 B disk is somewhat smaller than the 5.9 AU derived by 



Fur lan et al.1 (120071 ). Given the differences in the adopted stellar parameters, dust populations, and 
the prescribed physical structure of the cavity edge between these two studies, this apparent i? cav 
dis crepancy is p r obabl y irrelevant. However, if the cavity extended out to the radius determined 
by iFurlan et al.l (|2007l ) , we would expect to see a null in the 880 urn visi bilities on ~350-400 kA 
baselines, corresponding to angular scales of ~0 / . / 3 (see lHughes et al.l 120071 ). Since no such null is 
noted in the SMA data (see Figure [3]), a smaller cavity size is appropriate. 

The additional cavity parameters, {R m , <5s, 5h, /sg}, are not well-constrained and are some- 
what mutually degenerate. For example, one could adjust i?; n and compensate for the modified 
infrared emission by varying 5-% ■ Nevertheless, in practice the range of viable parameter-space is 
relatively small, and the values selected here are meant only to be illustrative and facilitate rea- 
sonable estimates of the disk edges. Both disks contain only a miniscule ~10 -10 M© of dust inside 
R ca v, a fraction / sg (by mass) of which is comprised of small grains. The scale-height modifications 
at -Rcav ($h) are important, but also modest, amounting to a change from 0.19 to 0.26 AU for 
HD 98800 B and 0.0026 to 0.0032 AU for Hen 3-600 A. The vertical structure parameters inferred 
for both disks are sim ilar to their non-truncated counterparts in the Ophiuchus star-forming region 
((Andrews et al.l l2009). but their total dust masses are an order of magnitude smaller than the typi- 
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cal T Tauri disk (I Andrews fc Williams!l2007bl ), Note that the Hen 3-600 A disk has a slightly higher 
Md us t than HD 98800 B, despite the latter being brighter at 880 /mi. The apparent discrepancy is 
a result of the significantly cooler dust temperatures produced around the later-type Hen 3-600 A 
stars. 

Because there is no spatially resolved information in the SEDs alone, those data can be repro- 
duced equally well with essentially any value for the outer disk edge. Not surprisingly then, the SED 
models for different -R ou t values shown in Figure [2] are practically identical. However, the synthetic 
880 fim continuum visibilities computed from the same disk models exhibit distinctive signatures 
for different outer edge locations. Figure [3] displays the azimuthally- averaged representations of 
the observed and synthetic visibilities for t he HD 98800 B and Hen 3-600 A disks, deprojected to 
account for the disk viewing geometries (see lLay et al.lll997l ). The real parts of the SMA visibilities 
decrease with deprojected baseline, demonstrating that both disks are resolved. Moreover, the 
imaginary parts of the visibilities are essentially zero, indicating that there is no strong evidence 
for large departures from axisymmetry. To facilitate a fair comparison of the data and models, 
the model visibilities are first convolved with a / /15 Gaussian to compensate for the atmospheric 
decoherence discussed in §2. After this "seeing" correction, it is clear in Figure [3] that the SMA 
data are best reproduced for i? out values of ~10-15AU for the HD 98800 B disk and ~15-25AU 
for the Hen 3-600 A disk. These values could be smaller if the decoh erence is under-estimated, b ut 
they could only be larger if S were a steeper function of radius (see lAndrews fc Williams! l2007al ) . 



5. Discussion 



Through radiative transfer modeling of resolved 880 /im continuum emission and pan-chromatic 
spectral energy distributions, we have characterized the key properties of the circumstellar material 
in two of the nearest pre-main sequence multiple star systems. In principle, the disk structures 
determined here can be used to help assess the potential for planet formation in the presence of a 
stellar companion and to test theoretical calculations of star-disk tidal interactions. The latter task 
has further implications for disk studies, as empirical constraints on these dynamical interactions 
can be used as high mass ratio (/i) touchstones for modeling the impacts of young giant planet 
companions on their local disk environments. 

Unf ortunately, very little is known about the stellar orbits in the Hen 3-600 hierarchical triple 
system. iMuzerolle et al.l (120001 ) noted that the primary appears to be a double-lined spectroscopic 
binary with a large radial velocity difference (~40km s —1 ), although repeated measurements have 
not been able to extract a reliable set of orbital parameters ([Torres et al.ll2003l ). If the inner trunca- 
tion radius estimated for the Hen 3-600 A disk, i? C av ~ 1 AU, i s indeed set by the dynamical clearing 
effects of the central spectroscopic binary, the calculations by I Artvmowicz Lubowl ()1994l ) predict 
that the orbital separation should be in the range a ~ 0.2-0.6 AU, depending on the eccentricity 
(mass ratio-related effects are small when /x ~ 0.1-0.5). The visual companion Hen 3-600 B has 
been located at angular separations of 1.4-l'.'5 to the southwest (PA ~ 215-230°) of the spectro- 
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scopi c binary system, corresponding to projected spa t ial separations of 63-67 AU (jde la Reza et al 
19891 : iReipurth & Zinneckerl Il993l ; IWebb et al1ll999l ; brandeker et alJbood ). While there are no 
available constraints on the A-B orbit, projection effects should only produce up to a ~20% uncer- 
tainty in the physical separation if the disk plane and orbit are aligned. Assuming that the similar 
spectral types of the A and B components imply a roughly equal-mass system (// ~ 0.5), tidal 
intera ctions would externally trun cate the Hen 3-600 A disk at a radius of ~25AU for a circular 
orbit (jArtymowicz &; Lubowlll994h . The outer disk edge inferred here, R out ~ 15-25 AU, is in good 
agreement with that prediction, with the uncertainty permitting moderate orbital eccentricities. 

The stellar orbital parameters in the HD 98800 quadruple system are comparatively more 
robust, thanks to long-term and intensive monitoring efforts. Most recently, iBoden et al.l (12 005) 



combi ned infrared interferometric observations with the extant spectroscopic data from 



Torres et al 



()1995l ) to measure the orbital elements and indivi dual stellar properties in the HD 98800 B spec- 
troscopic binary. According to the simulations of IPichardo et al.l ( 20081 ) . this nearly equal-mass 



(n = 0.43), eccentric (e = 0.78), and close (a = 1.05 AU) system can harbor a circumbinary 
disk with an innermost stable orbit at a radius of ~3.7AU. A similar value can be deduced 



fr om the calculations of iHolman &: Wiegertl (|1999i ). or from a slight extrapolation of the results 
of I Artvmowicz fc Lubowl (|1994i ). That value is practically identical to the internal truncation ra- 
dius determined from the inf r ared s pectrum morphology, R cav = 3.5 AU. For the wider HD 98800 
A-B visual pair, iTorres et al.l (|1995l ) compiled a relative astrometric record of the resolved system 
over the past century, along with a single-epoch measurement of the systemic radial v elocity dif- 
ferenc e between the components. With that information and the Hipparcos parallax, iTokovinin 
(|1999l ) determined a family of acceptable A-B orbits with ~300-430 yr periods, moderate eccentric- 
ities (e = 0.3-0.6), and a ~ 50-80 AU. For referenc e, the relative geometries of the HD 98800 B disk 
and an example A-B orbit from ITokovinin! ( 1999 ) are shown together in Figure EJ projected onto 
the plane of the sky. Note that the A-B orbital plane (i w 88°, $7 « 185°) is not aligned with the B 
spectroscopic binary and disk plane (i ~ 67°, $7 ~ 338°), making a quantitative prediction for the 
external disk truncation difficult. Although the aforementioned theoretical calculations all assume 
coplanar interactions, we can use them as a guide to assess whether or not the outer truncation ra- 
dius determined here is reasonable. Assuming that the A a nd B spe ctroscopic binar i es are roughly 
equal in total mass, the models of I Artvmowicz &: Lubowl (|1994l ) or IPichardo et al.l (|2005l ) predict 
that the HD 98800 B disk should be truncated at a radius of ~8-20 AU, depending on the assumed 



disk viscosity and which ITokovinin! (jl999i ) orbit is adopted. Those values are indeed consistent with 
the estimates based on modeling the resolved SMA data, where i? ou t ~ 10-15 AU. 

In general, the locations of the disk edges derived from our radiative transfer modeling are in 
good agreement with the theoretical predictions of disk truncation due to star-disk tidal interac- 
tions. However, there are some notable problems in the details that should be brought to attention. 
One of these involves the proper interpretation of the inner disk radius, R[ n , and the small amount 
of dust between it and the cavity edge (R cm )- F° r the disk models used here, that small dust mass 
is necessary to reproduce the weak infrared excesses and silicate emission features present in the 
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observed SEDs. However, the dynamical models suggest that there are no long-term stable dust 
orbits inside the inferred R CSuV . Perhaps the simplest way of reconciling this discrepancy is to argue 
that the adopted parameterization of the inner disk structure is inappropriate or incomplete. For 
example, similar infrared excess features could potentially be reproduced by modifying the dust 



conte nt, shape, or symmetry of the disk at the cavity edge (e.g.. Ilsella Sz Nattall2005l ; lAkeson et al 



20071 ) . With the same data for the HD 98800 B disk, iFurlan et al.1 1)20071 ) employed a narrow ring 
interior to their slightly larger cavity edge to reproduce the infrared spectrum with equivalent ac- 
curacy. The ability to obtain quality fits to the observations with such different interior structures 
highlights the model degeneracies, which will remain until the inner disk edge is spatially resolved. 

But if the inner disk structures determined here are taken at face value, the material in- 
ferred to be inside the cavity edge could be interpreted as dynamically unstable, perhaps repre- 
sentative of a low-dens it y inward flow due to Poynting-Robertson drag or disk accretion (e.g., see 
Pierens fc Nelson! l2008l ). IFurlan et al.1 (|2007l ) suggested the former mechanism to explain the nar- 
row ring from 1.5-2 AU they inferred in the HD 98800 B system, but rightly acknowledged that the 
empty region from 2AU to their R cav (5.9 AU) was incompatible with such a flow (note also that 
a significant fraction of such a ring would be inside the spectroscopic binary orbit). The structure 
parameterization adopted here is continuous from i? cav inwards to i?j n (which corresponds to the 
spectroscopic binary apastron distance), and so avoids this complication. Therefore, P-R drag 
could potentially account for the material inside the stable orbit radius in this disk if there is a 
continually replenished dust supply near i? cav and viscous gas drag forces are negligible. At least 
the latter condition appears to be me t for the HP 98800 B disk, where n o evidence for accretion 
or molecular gas has been uncovered (IDent et alJl2005l ; ISalyk et al.ll2009l ). The same is not true 



for the Hen 3-600 A disk, which exhibits the clear signatures o f a M ~ 5 x 10 11 yr 1 gaseous 



accretion flow from the outer disk to the spectroscopic binary (jMuzerolle et al.ll2000l ). In this case, 
the small amount of material inside R cscv may be tracing dust grains swept up in that flow. 

When future astrometric and spectroscopic monitoring has fully determined the visual orbits 
in the HD 98800 and Hen 3-600 systems, the lack of millimeter emission around the secondaries 
can also help constrain how their circumstellar material evolved. In principle, the tidal truncation 
calculations described above predict similar (pe rhaps slightly smaller) outer edge locations for th e 
disks around the nearly equal- mass secondaries (|Artymowicz Sz Lubowill994l : lPichardo et al.ll2005l ). 
If the individual disks initially had the same masses, the similar disk sizes would have produced com- 
parable millimeter emission levels from both components. Obviously this is not the case for these 
targe ts, nor for a handful of other resolved multiple systems ( Jensen Sz Akeson 2003 1; Patience et al. 
2008). This could be an artifact of the multiple star formation process itself, as iBate Sz Bonnel] 



(|1997l ) have argued that compe titive accretion will produce higher initial disk masses around pri- 
mary stars (see also iBatd |2000j) ■ Alternatively, the dynamical interactions in these systems could 
be considerably more complicated than the standard tidal truncation calculations have assumed, 
with the secondary disks suffering more severe disruption. For example, the interactions in non- 
coplanar systems like the HD 98800 A-B visual binary could be catastrophic, producing a much 



more complex circumstellar structure (jVerrier &c Evansll2008l ) 



The Hen 3-600 and HD 98800 multiple star systems should continue to be exploited as essential 
test cases for evaluating how external dynamical interactions can contribute to disk evolution, and 
therefore affect the potential for planet formation in the presence of stellar companions. The 
unique proximity of these particular systems is crucial, providing both high sensitivity and spatial 
resolution to circumstellar dust structures as well as rare opportunities to determine precise orbits 
for stellar pairs with different separations within the same system. While this study has tackled the 
initial challenge of resolving the millimeter emission in these systems and relating it to the global 
disk structures, the real utility of these targets now lies with continued work that focuses on the 
details. For example, future observations of these disks with the Atacama Large Millimeter Array 
could be sensitive to the structu ral asymmetries or recently ejected materia l streams expected from 



non-coplanar tidal interactions (jAkeson et al. 



2007 



Verrier fe Evans! 120081 ). However, it is equally 



important to calculate the stellar orbits in these nearby systems, particularly in the less-studied 
case of Hen 3-600. In fact, there are exciting opportunities to simultaneously refine the stellar orbits 
and directly probe star-disk interactio ns for the HP 9 8800 system, as the visual A-B binary will 
reach periastron in the next ~15 years (jTokovininl 1 1 999l ) . Therefore, with continued monitoring the 
circumstellar material in the Hen 3-600 and HD 98800 systems can eventually serve as the primary 
observational calibrators for models of star-disk (as well as planet-disk) tidal interactions. 
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Aa ["] Aa ["] 

Fig. 1. — Aperture synthesis images of the 880 (340 GHz) dust continuum emission in the 
HD 98800 (left) and Hen 3-600 (right) multiple star systems. Contours start at 3 a (lOmJy beam -1 ) 
and are shown at 2 a intervals. The synthesized beam dimensions are marked in the lower left corner 
of each image, and a 30 AU scale bar is included in the upper left. The images are centered at the 
peak of the continuum emission, corresponding to the positions of the HD 98800 B and Hen 3-600 A 
spectroscopic binaries. Each stellar component is marked with a blue star symbol. 
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Fig. 2. — Broadband spectral energy distributions for the HD 98800 B (left) and Hen 3-600 A (right) 
disks, constructed from measurements in the literature. Light gray curves mark the composite 
stellar photospheres adopted for the radiative transfer modeling, while the black curves show the 
Spitzer IRS spectra corrected for the contamination of the nearby companions (see text). The 
overlaid curves are disk structure models for different values for the outer disk radius, i? utj as 
labeled in the lower left corner of each plot. The individual models are difficult to distinguish here, 
but show clear differences in terms of the 880 /jm visibilities (see Fig. [3]) . 
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Fig. 3. — The azimuthally-averaged 880 continuum visibilities (real and imaginary parts) as a 
function of the deprojected baseline length for the HD 98800 B (left) and Hen 3-600 A (right) disks. 
The overlaid curves are truncated disk structure models with different values for the outer edge 
radius, -R ut> as labeled in the lower left corner of each plot. 
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Fig. 4. — The inferred geometry of the HP 98800 B disk (grayscale) is compared with the A-B visual 
binary orbit (red) computed by iTokovininl (|1999l ). projected onto the plane of the sky. The two- 
toned gray shading of the disk corresponds to the quoted range of possible R out values (~10-15 AU). 
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Table 1. Adopted Stellar Parameters 



Star 




R* 


log 5 


Ay 




[K] 


[R©] 


[cm s -2 ] 


[mags] 


HD 98800 Aa+Ab 


4500 


1.30 


4.00 


0.0 


HD 98800 Ba 


4250 


1.16 


4.25 


0.3 


HD 98800 Bb 


4000 


0.90 


4.25 


0.3 


Hen 3-600 Aa+Ab 


3350 


1.20 


4.50 


0.0 


Hen 3-600 B 


3350 


1.00 


4.50 


0.0 



Note. — The effective temperatures, radii, gravities, 
and extinctions determined for each stellar component, 
based on matchin g resolved optical /nea r- infrared pho- 
tometry to scaled iLejeune et al.l (119971 ) spectral tem- 
pla tes. The parameter s for HD 98800 B are based on 
the lBoden et al.1 (120051 ) measurements: the templates 
were summed into a composite spectrum for the radia- 
tive tr ansfer modelin g. The Ay estimates correspond 
to the iMathisI (|l990i ) extinction law. 
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Table 2. Estimated Disk Parameters 



Parameter 


HD 98800 B 


Hen 3-600 A 


Rin [AU] 


2.0 


0.4 


Rc. v [AU] 


3.5 


1.0 


Rout [AU] 


10-15 


15-25 


M dust [M ] 


3 x 10~ 6 


7 x 10~ 6 




~1000 


~500 


H W av [AU] 


0.5 


0.3 




0.10 


0.07 




1.4 


1.2 


/sg 


0.25 


0.10 


i [°] 


67 


36 


PA [°] 


158 


169 



Note. — Disk structure parameters esti- 
mated from the radiative transfer modeling. 
The parameter definitions are given in §4. 
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